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ABSTRACT 

The origin of the diffuse extragalactic gamma-ray background (EGB) has 
been intensively studied but remains unsettled. Current popular source candi- 
dates include unresolved star-forming galaxies, starburst galaxies, and blazars. 
In this paper we calculate the EGB contribution from the interactions of cos- 
mic rays accelerated by Type la supernovae (SNe), extending earlier work which 
only included core-collapse SNe. We consider Type la events not only in star- 
forming galaxies, but also in quiescent galaxies that lack star formation. In the 
case of star-forming galaxies, consistently including Type la events makes little 
change to the star-forming EGB prediction, so long as both SN types have the 
same cosmic-ray acceleration efficiencies in star-forming galaxies. Thus, our up- 
dated EGB estimate continues to show that star-forming galaxies can represent 
a substantial portion of the signal measured by Fermi. In the case of quiescent 
galaxies, conversely, we find a wide range of possibilities for the EGB contribu- 
tion. The dominant uncertainty we investigated comes from the mass in hot gas 
in these objects, which provides targets for cosmic rays; total gas masses are as 
yet poorly known, particularly at larger radii. Additionally, the EGB estimation 
is very sensitive to the cosmic-ray acceleration efficiency and confinement, espe- 
cially in quiescent galaxies. In the most optimistic allowed scenarios, quiescent 
galaxies can be an important source of the EGB. In this case, star-forming galax- 
ies and quiescent galaxies together will dominate the EGB and leave little room 
for other contributions. If other sources, such as blazars, are found to have impor- 
tant contributions to the EGB, then either the gas mass or cosmic-ray content of 
quiescent galaxies must be significantly lower than in their star-forming counter- 
parts. In any case, improved Fermi EGB measurements will provide important 
constraints on hot gas and cosmic rays in quiescent galaxies. 
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Introduction 



The first observation of th e diffuse extragalactic ga mma-ray background (EGB) was 
reported by the SAS-2 satellite (IFichtel et al.lll977l Il978[ ). Recently, the Fermi Gamma-ray 
Space Telescope up_dated the EGB de termination from the Energetic Gamma-ray Experi- 
me nt Telescope ([Sreekumar et al. 1998 ) and provided the most reliable EGB observations so 
far ( Abdo et al. l2009al ). Measurements of EGB are model depen dent in that they r equire 
subtraction of the large foreground emission from our Galaxy (e.g.. iHunter et al.lll997l ). The 
accuracy of the EGB measurement thus greatly depends on our understanding of the Galactic 
emission. Despite the difficulty in its observation, the EGB encodes important information 
about the highest-energy environments in the cosmos. 

The EG B arises from t he combination of all th e unresolved extragalactic gamma-ray 
sources (e.g.. lDermerll2007bl ; IStecker k Ventersll201ll ). "Guaranteed" EGB components arise 



from unresolved counterparts of known extragalactic populations, nam ely, blazars (those ac- 



tive galactic nuclei t hat have their relativistic je t s pointing at us, e.g. . iPadovani et al 



Stcckcr et al. 



Venters 



1993 



Muk heriee k Chiang 



2010; Venters k Pavlidou 



2011 



galaxies and starburst galaxies (e. 



1999 



Pavlidou k Venters 20081 



Inoue k Totani 



Dermer 



1993 



2007a 



200 



Pavlidou k Field 



20061 : iThompson et al.ll2007l : ISteckerll2007l : Fields et al.ll2010l : lMakiva et al.ll2011l ). Additional 



si), a s well as normal star-forming 



2001 



2002 



Prodanovic k Fields 



EGB contributions might arise from more exot ic sources, such as dark matter annihilation 
( ISilk k Srednickil Il984j ; iRudaz k Steckerl Il991 ) ; anni hilations at the boundaries of cosmic 



matter a nd antimatter domains (IStecker et al.l Il97ll ) , massive black hol es at redshifts of 
z ~ 1 00 (IGnedin k Ostrikerlll992f ). and primordial black hole evaporation (IPage k Hawking 



19761 ). 



In this paper we will focus on the EGB contribution from both star-forming galaxies and 
quiescent galaxies. Quiescent galaxies refer to galaxies with little or no active star formation, 
and these objects have not been included in EGB estimations. In terms of galaxy types, 
quiescent galaxies usually include all elliptical galaxies and some SO galaxies. However, the 
important factor for the EGB estimation is not the galaxy type but the amount of star 
formation. Therefore we will separately consider star-forming and quiescent galaxies, and 
assume no star formation in quiescent galaxies. We will not consider the EGB contribution 
from starburst galaxies in this paper, due to the larger uncertainty in t he cosmic-ray prop- 



agation in such galaxies (e.g., IThompson et al.l 120071 ; lLacki et al.l 1201 ll ). Also, recent work 



sugge sts that starburst galaxies only have small contribution to the EGB (IStecker k Venters 
20111 ). We follow the criterion adopted in lFields et al.l (120101 ) to distinguish star- forming and 
starburst galaxies. The EGB energy range we consider in this paper is from ~ 30 MeV to 
~ 30 GeV, which is the energy range covered by Fermi data and includes the regime in which 
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star- forming galaxies may contribute substantially to the EGB. 



Our focus here is on the EGB contribution arising from hadronic cosmic-ray interac- 
tions with the interstellar medium (ISM) of their host galaxies, specifically pion production 
and decay pp — > tt° — > 77. The most favored possibility of the cosmic-ray production 
sites in galaxies is supernovae (SNe). Fermi and ai r Cerenkov observations detect indi 



vidual remnants of both core-collapse ( hereafter CC: lAbdo et al. 



Weekes et al.lll989l ) and Type la events (lAcciari et al.ll2011 



2010 



Acero et al. 



Acciari et al. 2010 



20101 ). The energet- 



ics of these objects are consistent with the require ments of efficient cosmic-ray acceleration 



(e.g., 



Abdo et al.l l2010t iReynolds &: Ellison! Il992l ) , and the Ge V spect ra of these objects 



Abdo et al. 2009b 



are consistent with pionic emission an d thus hadronic acceleration (e.g. 
Tanaka et al .1120 111 ; Ellison et al.l 120121). Diffuse Galactic emission is similarly consistent with 
pionic emission dominating (jAbdo et al. 2009al Jcl). These data thus give empirical grounding 



to the long-held b elief that SNe of all typ e s are the dominant engines that accelerate Galacti c 



Baade fc Zwickvlll934j ; iGinzburg & Svrovatskiilll964j : Ellison et al.lll997h . 



cosmic rays (e.g. 

We are thus interested in SNe of both types, all of which accelerate cosmic rays in their host 
galaxies. 

Many groups have st udied the EGB emission from cosmic rays accelerated by SNe in 



star-forming galaxies (e.g..lDar fc Shayivlll995l ; iProdanovic fc Fieldsll2006t Fields et al.ll2010 



Stecker fe Venters! 1201 It iMakiya et al.ll201ll ). Som e estimations sugg est that star-forming 
galaxies can be the dominant source of the EGB ( [Fields et al.l 120101). while other groups 



predict that a major contribution of the EGB comes from blazars (jStecker & Venters! 12011 



Makiya et al.ll201ll ; llnoue fc Totanill2009l ). However, there exist large uncertainties from the 
source inputs. Most of the analyses regarding star-forming galaxies focused on the EGB 
contribution from cosmic rays accelerated by CC SNe and implicitly assume that only these 
ev ents accelera t e cosm ic rays. We extend the analysis of the EGB from star-forming galaxies 



m iFields et al.l ( 120101 ) to include Type la SNe as accelerators in the Milky Way and in other 
galaxies. 

CC SNe arise in massive stars with short lifetimes, and thus trace ongoing star formation. 
In contrast, Type la SNe result from thermonuclear runaway of white dwarfs accreting mass 
from their companion stars and hence are related to star formation with some delay time. For 
this reason, observations have shown that Type la SNe exist in both star -forming galaxies 



and quiescent galaxie s, while CC SNe are rarely seen in quiescent galaxies (jFilippenkol 12001 



Mannucci et al.ll2005l ). Observations have suggested that the intrinsic cosmic C C SN rate is 



abou t 5 times higher than the intrinsic cosmic la SN rate at redshift z < 0.4 ( jBazin et al. 



20091 ). Also, studies suggest that the la rate in a star-forming galaxy is much larger than 
that in a quiescent galaxy; this reflects the distribution of delay times between progeni- 
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tor birth and la ex plosion, which is weighted towards short delays ( IMannucci et al.l 12005 
Sullivan et aUbood l. 



The efficiency of cosmic-ray acceleration by SNe remains poorly understood but is crucial 
for understanding cosmic-ray acceleration physics as well as SN energy feedback. Theories 
propose that cosmic rays are produced by diffusive shock accelerat ion in the blast waves from 
SN explosions (e.g., ISchlickeiserl Il989l : iBerezhko fc Ellison] Il999h . Current studies suggest 
that ~ 30% of the initial kinetic energy from a SN needs to be transferred to cosmic-ray 
acceleration if we assume that SNe are the dominate sources for cosmic- ray production an d 
the nucleosynthesis of lithium, beryllium, and boron in the Milky Way ([Fields et al.ll200ll ). 
Also, some theoretical predictions expect the cosmic -ray acceleration effic iency in quiescent 
galaxies is much lower than in star-forming galaxies. iDorfi Voelkl ( 119961 ) suggest that only 
< 1% of the total explosion energy goes into cosmic-ray energy in quiescent galaxies. This 
is because a SN blast will have weaker shocks due to the large sound speeds of the hot, 
low-density ISM in an elliptical galaxy. 

Understanding the SN rate and their efficiency in producing cosmic rays is critical for 
studying the EGB contributions from these galaxies. Our observational understanding of 
cosmic SNe will increase significantly when the next generation optical survey telescope, 
the Large Synoptic Survey Telescope (LSST), comes online during the next decade. LSST 
is planning to sca n the whole availa ble sky, repeated every ~ 3 days, with unprecedented 
survey sensitivity (llvezic et al.l 120081) . The project will observe ~ 10 5 CC SNe per year out 
to redshift z ~ 1 (ILien &: Fields! 120091 ) and ~ 5 x 10 4 Type la events out to redshift z ~ 0.8 
faailev et al.ll2009h 7 The cosmic SN rate in different galaxy classes can thus be measured 
via direct counting to high redshift with extremely low statistical uncertainty. 

In this paper, we will first describe the general formalism of estimating the EGB from 
cosmic rays accelerated by SNe in both star-forming and quiescent galaxies (§ [2D - We will 
then discuss the cosmic Type la rate in each galaxy classification that will be used in our 
EGB analysis (§[3]). The estimations of the EGB contribution from Type la SNe in star- 
forming and quiescent galaxies are presented in § H] and § 0, respectively. Additionally, we 
discuss the uncertainties in the EGB predictions in § Finally, we summarize the results 
in §0 



2. General Formalism 



The formalism we adopt generalizes that of IFields et al.l (120101 ) to account for both 
SN types. Integration of the gamma-ray contributions from each unresolved extragalactic 
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source over the line of sight to the cosmic horizon gives the well-known express for the EGB 
intensity, 



dl c 
dE ~ 4^r 



£ 7 (£ e m, z) (1 + z) 



dt 



dz, (1) 



dz 

where C 7 (E em , z) = dN ggbmraaj / dV cora dt em dE em is the comoving luminosity density (or emissiv- 



l-i 



ity) at rest-frame energy E em , and \dt/dz\ = [(l+z)H(z)) 1 = [{l+z)H ^JVL m (l + z) 3 + f^y] 
for the standard ACDM cosmology We use VL m = 0.27, fi A = 0.73, and H Q = 70 km s" 1 Mpc" 1 



from the seven-year Wilkinson Microwave Anisotropy Probe ( WMAP) data (IKomatsu et al. 



201lh . 



Because the pionic gamma-ray emission is produced from the interaction between cosmic 
rays and the hydrogen atoms in the ISM of each galaxy, the luminosity density is given by 
the product 

Cry = (L-y H. ga l), (2) 

of the pionic gamma-ray luminosity L 7 of an individual galaxy, times the galaxy number 
density, appropriately averaged. Our problem then divides into two parts. First, we must 
express a galaxy's gamma-ray luminosity L 7 in terms of galaxy properties such as SN rate 
and gas content, and relate these to galaxy observables. Then we must construct a luminosity 
function dn ga \/dL 1 for gamma-ray-emitting galaxies. 

We first turn to the pionic gamma-ray luminosity from an individual galaxy. This can 
be written as 



Lj(E em ) — J T 7T o^ 11 (E em ) nu dVi S M (3) 
= r>^ 77 (£ cm ) Ah (4) 

where T 7r o_ > . 77 (i? em ) represents a spatial average of the gamma-ray production rate per inter- 
stellar hydrogen atom. The total number jV H = J n H dVisM of hydrogen atoms in a galaxy 
is obtained by integrating the number density of hydrogen atom nu over the ISM volume. 
A/"h is proportional to the total gas mass M gas in a galaxy and can therefore be expressed as 
Ah = X-r M gas /m p , where Xn is the mass fraction of hydrogen atoms and m p is the proton 
mass. 

We take SNe (of both types) to be the engines of cosmic-ray acceleration; this implies 
that the cosmic-ray flux scales as $ cr oc A esc i?sN,eff • Here i?sN,efi is an effective SN rate 
weighted by the cosmic-ray acceleration efficiency e, discussed below. A esc is the escape path 
length, which quantifies the cosmic-ray confinement in a galaxy. In this paper we assume 
A esc to be universal and constant, which leads to a universal galactic cosmic-ray spectrum 
that is the same as that of the Milky Way. Thus the pionic gamma-ray production rate per 
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hydrogen atom in a galaxy should scale (E cm ) oc $ cr oc A esc -RsN,eff- We normalize 

the cosmic-ray spectrum to a known galaxy, which would be the Milky Way in our case, 
finding 

-RsN.efT 



rV — Ky"y (-^e 
yj-0 — y-yy V 



^cr 
<J>MW 



pMW • 

- rL SN,efT 



The pionic gamma-ray luminosity of a particular galaxy is thus 

\ _ pMW / rp \ 



pMW Ah ~ • 

-"-SN.efif '"P 



(5) 



(6) 



Because of their short lifespans, the rate of CC SNe (and short-delay la events) traces that 
of star formation: i?sN oc if). Thus we expect a star- forming galaxy's gamma-ray luminosity 
to scale M gas R SN oc M gas -0. The new class of Fermi-detected star-forming galaxies 



is consistent with this trend ( iLenain fc Walter! 1201 lh . 



This pionic gamma-ray spectrum always has a peak at E^m = at whi ch the two 

gamma-ray photons inherit the rest-mass energy of the decayed ir° (jSteckerl Il97ll ) . At large 
energy, the spectrum shows the same asymptotic index as that of the cosmic-ray spectrum, 
which we take to be 2.75. 

Both CC and Type la events should produce cosmic rays and hence pionic gamma 
rays. Therefore, the effective SN rate i?sN, e ff in Eq. (EJ) is a combination of the effective 
Type la rate R^s = £i a Ria. an d the effective CC SN rate -Rcc.cfr = ecc -Rcc; where 
ei a and ecc are the cosmic-ray production efficiencies of Type la and CC SNe, respectively. 
There exist different definitions of the cosmic-ray acceleration efficiency in current literature. 
For example, some studies present the efficiency as the fraction of the to tal cosmic-ray 
production energy out of the total kineti c energy output from a SN (e.g., iDorfi fc Voelk 
19961 ; Fields et al.ll200ll ; iHelder et al.ll2010l ). while other studies define the parameter as the 
perce ntage of the energ y flux that becomes relativistic particles after crossing the shock 
(e.g., Ellison et al.ll2007l ). Most of these definitions describe the fraction of the SN explosion 



energy transferred to cosmic rays. Here, we define the cosmic-ray acceleration efficiency e 
to be the ratio of the SN baryonic explosion energy £"sn to the resulting cosmic-ray energy 
E cr , i.e., e = E^/Es^. Therefore, if we assume that all SNe have the same explosion energy 
and the produced cosmic rays have the same energy spectrum, the cosmic-ray acceleration 
efficiency will be proportional to the total cosmic-ray production in a galaxy over the SN 
rate in that galaxy, i.e., e oc $ cr /(A esc -Rsn) 



For the Milky Way, then, we have -RsN^, ff 



dMW i , pMW 
eia,MW -Kla + e CC,MW -^CC ■ 



Since we normalized our prediction to the gamma-ray production in the Milky Way 
(Eq. [6]), the important factor in the calculation is not the absolute value of e, but the 
difference between the acceleration efficiency e in different SN types (la and CC) and galaxy 
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classes (quiescent and star- forming). Specifically, we will need to specify the ratios ei a /ecc 
and 6Q/es- Unfortunately, these two fractions are poorly known. Thus for our fiducial 
numerical results, we will take eiaAcc = 1 an d cq/cs = 1- Furthermore, we are unaware 
of any evidence for a substantial difference between the cosmic-ray acceleration efficiencies 
between la and CC SNe, and thus we will hereafter drop the la and CC notation in the 
acceleration efficiency to simplify the discussion. However, we will retain the notations of 
the acceleration efficiencies for different galaxy type eq and es in our formalism as an explicit 
reminder that the efficienc ies are likely to depend on galaxy environ ment, as expected by 
some theoretical analyses (IDorfi &: Voelklll996r iHein fc Spanierl 120081 ) . Further possibilities 
of choosing different cosmic-ray acceleration efficiencies will be discussed in § 16.21 

Star-forming galaxies contain both Type la and CC SNe. Their pionic gamma-ray 
luminosity density £ 7i s can be calculated by averaging over the galaxy density n galaxy , 

eS-Rla,S ^galaxy) + (M gas €sRqC ^galaxy) 



-7,S 



■pMW 

J" ^0 V 



77 



(zp \ H (M gas 



(7) 



In quiescent galaxies, there are almost no star formation. We will assume the star-formation 
rate (and thus the CC SN rate) to be zero in a quiescent galaxy. However, Type la SNe do 
exist in quiescent galaxies because these events occur some time after the star formation. 
Therefore the pionic gamma-ray luminosity density in quiescent galaxies £ 7i q only comes 
from Type la events, 



F 7r o_ i , 77 (-E'em) Xn 



c 



7,Q 



e Q0 R?a 

Fjj-O^yy ("^e 

standard quiescent galaxy Q0 for normalization. However, since no gamma-ray emission 
from a (jet-less) quiescent galaxy has ever been measured, we will still adopt the values of 
the Milky Way and estimate the gamma-ray luminosity density for quiescent galaxies as 

Xn (M g as e Q ^Ia,Q ^galaxy) 



rrir. 



(M g 

as Rla.,0, ^galaxy)- 



(8) 



and are the gamma-ray production rate and Type la event rate in a 



7,Q 



r 



MW 



>77 



™ c pMW I , dMW • 



(9) 



Note that since the gamma-ray production from the Milky Way comes from both Type la 
and CC SNe, r^ w v . (E em ) needs to be normalized to the total SN rate in the Milky Way 



instead of just the Type la rate. 

The total pionic gamma-ray luminosity density will be a combination of emissions from 



7, tot 



The EGB contribution 



both star-forming and quiescent galaxies, that is, C 
from cosmic rays accelerat ed by CC SNe has been c arefully examined i n lFields et al.l (120101 ) 



and by other groups (e.g., IStecker & Venters! 1201 ll ; iMakiya et al.ll201ll ). Here we will focus 



on the EGB contributions related to Type la events. In our calculatio n, we do not include 
the intergalactic EGB absorption > 30 GeV fjSalamon &: Steckerlll998l ). 
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The Cosmic Type la Supernova Rate in Star-forming and Quiescent 

Galaxies 



Type la SNe do not all trace ongoing star formation, because these events have different 
origins from CC SNe. The prevailin g scenarios for Type la SN origin include merging of 
two white dwarfs (double degenerate. IWebbinklll984l ). o r a white dwarf accreting from mass- 
overfl ow of its supergiant companion (single degenerate, iNomoto et al.lll984[ llben fc Tutukov 



19841 ) . Both of these scenarios involve white dwarfs merging in a binary system, and thus 
Type la SNe are delayed from the formation of the progenitor stars. For this reason, Type 
la SNe are found in all galaxies, including the quiescent galaxies where there is no longer 
star-forming activity. A complete account of the Type la SN contribution to the EGB must 
therefore include contributions from events in star-forming and quiescent galaxies. 

There are many studies of the comoving cosmic la rate density TZi a = dNi a / (dV com dt) 



as a function of redshi 



Scannapieco &: Bildsten 



t, most of which 



2005; 



Sullivan et al 



bcus on the distribution of delay times (e.g., 



20061 : Kuznetsova et al.ll2008l : iDildav et aLBOlOb 



Horiuchi fc Beacoml |2010| ; iGraur et al.ll201ll ). Our adopted total (all-galaxy) cosmic Type 



la rate is based on current observational data. To find this we use the best-fit redshift and 



delay-time distribution of iHoriuchi fc Beacoml (120101 ). integrated over all delay times (black 
curve in Fig. [T]). Although the number of observed cosmic la events is rapidly increasing, 
the data remain sparse beyond z ~ 1, where the rates are thus poorly constrained. 

The cosmic la rate in different galaxy classes (i.e., star- forming and quiescent galaxies) 
as a function of redshift is still poorly understood. In our calculation, we adopt a con- 
stant value for the cosmic la rate in quiescent galaxies out to redshift z ~ 2 (red curve 
in Fig. [T]). We use a normali zation based o n obse rvations of la rate per stellar mass in 



quiescent galaxies provided by Sullivan et al.l (120061 ) and a non-evolving stellar-mass func- 
tion from iPannella et al.l (120091 ) (see detailed discussion in § EJ). Recently, Type la SNe 



within z < 1 ( 


Gal- Yam et al. 


2002 


Sharon et al. 


2007 


2008; 


Dildav et al. 


2010al: 


Sharon et al. 


2010 


). Since 



of quiescent galaxies, these results are consistent with our assumption of a constant la 
rate in quiescent galaxies. Moreover, these measurements show that la rate in clusters is 
~ 10~ 13 yr _1 Mq 1 , which is very similar to the value we adopted from I Sullivan et al.l ( 2006 ) 
0. However, while we expect a roughly constant cosmic la rate in quiescent galaxies for 
moderate redshift, this trend must fail at some redshift. To be conservative, we therefore 



1 We are greatly thankful for the anonymous referee for pointing out the measurements of la rate in galaxy 
clusters. 
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placed an artificial cutoff of the cosmic Type la rate at z — 2, beyond which uncertainties 
in the cosmic Type la rate observations remain substantial. 

The cosmic la rate in star-forming galaxies (blue curve in Fig. [TJ) can thus be obtained 
by subtracting the cosmic Type la rate in quiescent galaxies (red curve in Fig. [1]) from the 
total cosmic la rate (black curve in Fig. [1]). 

Figure [U shows the adopted cosmic la SN rate as a function of redshift. Although the 
uncertainty in the rate increases significantly at higher redshift, most of th e EGB from Type 



la SN e, like that from CC SNe, arises from events at lower redshift (z < 1) ( lAndo fe Pavlidou 



20091 ). For example, in our calculation ~ 50% (~ 70%) of the EGB flux comes from sources 
within z < 1 (z < 1.3). Therefore the choice of the Type la rate at z > 1 only has a small 
effect on the final estimation of the EGB. The solid black curve plots the total cosmic la 
SN rate in both star-forming and quiescent galaxies. The red curve shows the cosmic la 
SN rate in only quiescent galaxies. The blue curve represents the cosmic la SN rate in only 
star-forming galaxies. 

Figure [T] also shows the cosmic CC SN rate (dotted black curve), which is higher than 
the la rate by a factor ~ 5 at z ~ and increase to a factor of ~ 10 at z ~ 1. This 
immediately suggests that we should expect CC events to dominate the star-forming EGB 
signal, with the la contribution at a < 20% level. We will see that this is roughly the case 
for the la contribution from star-forming galaxies, but for la events in quiescent galaxies the 
situation is much more uncertain. 



4. The Extragalactic Gamma-ray Background from Type la Supernovae in 

Star-forming Galaxies 

As described in § |2j the EGB luminosity density £^ a s is dominated by two physics 
inputs: the SN rate in a galaxy, which is associated with the amount of cosmic rays, and the 
total gas mass of that galaxy, which accounts for the total hydrogen targets that interact 
wi th the cosmic rays . To reflect these two physics inputs, we follow the approach adopted 



m 



Fields et al.l (120101 ) and rewrite the EGB contribution from Type la events (the first term 



in Eq. [7]) as below, 



^7,S - " nMW , £ pMW 1 I rM 77 l fi cmJ \iW gaSjS /, (±U) 

fMW -n-ia + e MW -^CC m P 
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Fig. 1. — Adopted cosmic SN rate. The solid black curve plots the total cosmic Type la 
rate; blue curve plots the cosmic la rate in star-forming galaxies; red curve plots the cosmic 
la rate in quiescent galaxies. The cosmic CC SN rate is plotted as dotted-black curve for 
comparison. 



where 



/vf \ _ (^gas,S #Ia,S ™gala X y,s) ,,,, 
\-Kla,S n galaxy,S/ 



dn 



J dLn a:Z Ri^siLua, z) 



dn 



dLuc 



and 7£ia,s = (-Ria.s %ai axy ,s) is the cosmic Type la rate in star- forming galaxies, as shown in 
Fig. [TJ In a star-forming galaxy, the galaxy gas mass M gas g and the galaxy Type la rate 
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i?i a can be related to the star-formation rate in that galaxy, which can be connected to the 
observab le Ha luminos ity L Hq , z of the galaxy by ip(Ln a , z)/(l M Q yr _1 ) = L HQ)Z /(1.26 x 
10 34 W) (lHopkinsll2004l ). Therefore we express the gas mass M gas S and the Type la rate i?i a 
in terms of L^ a z . The correspondin g galaxy luminosity fun ction at this wavelength can be 
expressed by the Schechter function (INakamura et al.ll2004l ). 



At a specific redshift, the gas mass in star- forming galaxies M gas s and the star-formation 
rate can be connected by 



M gaSiS = 2.8 x 10 9 M (1 + z) 



1 M yr- 1 



(13) 



with (5 = 0.571 and u = 0.714, as shown in lFields et al.l f 20101 ) . The Type la rate in a galaxy 
can be linked to the star-formation rate via some delay-time distribution A(t), 



Ruiz) oc 



t(z) 



ip(t - r) A(r) dr, 



(14) 



where t(z) is the corresponding cosmic age at redshift z. The delay-time distribution A(r) 
gives the probability that a Type la SN explodes a time r after the progenitor's birth. More 
detailed discussion about the delay-time distribution can be found in Appendix [A] The 
galaxy luminosity function at a certain redshift for star-forming galaxies in the Ha band can 
be presented in the form of a Schechter function of 

dn n+, ( Lhq- 



dL 



H«,z 



(15) 



with a = 1.43 (INakamura et al.l 120041 ) . 



Because Type la SNe are delayed relative to star formation, the Type la rate in a 
galaxy depends on the star formation history of the galaxy via eq. (114p . Unfortunately, 
this past star formation rate is not easily determined for distant galaxies even when the Ha 
luminosity Lh«,z is available to give the instantaneous star formation rate. However, we can 
investigate the evolution in two simplified cases: pure luminosity evolution and pure density 
evolution. Pure luminosity evolution assumes that galaxy luminosities evolve with redshift, 
while galaxy density stays unchanged, i.e., in Eq. f[T5"j) has redshift dependence and n* 2 
does not. Pure density evolution assumes that galaxy density evolves with redshift, while 
galaxy luminosity remains constant, i.e., n* )2 in Eq. (|T5|) depends on redshift and L* iZ does 
not. The real situation should be bracketed by these possibilities. 

In either limit, we can take advantage of the fact that the well- measured cosmic star 



formation rate is /?* = (i/m ga \) ~ (-^Ha^-gai) ~ 
product L+zUi, 



K n* jZ . This fixes the redshift behavior of the 
so that in the limit where one factor is constant, the other must have the 



redshift dependence of the cosmic star formation rate. 
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4.1. Pure Luminosity Evolution 



In the case of pure luminosity evolution, there is no evolution of the galaxy density Thus 
in the Schechter function of eq. (fT5j) . n* jZ = ra* j0 , and thus all redshift dependence lies in Lh q , z . 
Therefore, evolution of the star-formation rate in each galaxy, and hence the evolution of the 
galaxy Ha luminosity -Lh«, z , must trace the general evolution of the cosmic star-formation 
rate p±. Under this assumption, we can show that (M gaSj s) is independent of the choice of the 
delay-time function (see derivation in Appendix |A]). When adopting the Schechter function 
for d " B r alaxy ' s , one will find that (M gas , s ) oc (1 + z)~ p oc (1 + z)~ p (^LV, with a 



local value of (M gaSj s)z=o = 6.8 x 10 9 M Q (see Appendix IA1 . 



>p*(*=o). 



The predicted EGB from Type la SNe in star-forming galaxies is plotted as the solid 
blue line in the left panel of Fig. |2j For comparison, the dashed blue line shows the EGB 
contribution from CC SNe in star-form ing galaxies. The shapes of the dashed blue lines 
trace the results in iFields et al.l (120101 ) . However, the normalization of the CC SN curves 
is lower by the fraction of the CC S N rate over the tota l SN rate (~ 0.8 from Bazin et al., 
2009), which is due to the fact that IFields et al.l (j2010[ ) have implicitly assumed that CC 
SNe produce all of the gamma-ray emission in galaxies. 

Figure [2] shows that the EGB from Type la SNe is around an order of magnitude lower 
than those from CC SNe, which is due to the lower Type la rate in star-forming galaxies. 
As noted above, this is easily understood as a reflection of the small Ia/CC ratio (Fig. [[]). 
The black curve in Fig. [2] presents the total EGB emission from both Type la and CC SNe 
in star-forming galaxies. Note that the to tal EGB emi s sion f rom star-forming galaxies turns 
out to be very similar to the prediction in IFields et al.l ( 2010 ). in which the authors assumed 
that all of the EGB contribution comes from the CC events. The reason is that even though 
we added the EGB con tribution from Typ e la SNe, we also lower the EGB emission from 
CC events estimated in Fields et al.l (120101 ) by the corresponding CC SN fraction. Also, the 
la to CC fraction does not change much within z ~ 1, which is the redshift range where most 
of the EGB signals originate. Hence the EGB contribution from CC SNe is always higher 
than those from la SNe by a similar factor. 

The shape of the EGB curves in Fig. [2] traces the general features of the pionic gamma- 
ray energy spectrum. This is because the observed EGB intensity at a specific energy 
originated from a combination of sources at different redshifts, as described in Eq. ([jrj. 
Therefore, the redshift evolution of the unresolved sources is smeared out in the energy plot 
and mostly affects the normalization of the EGB intensity but not the spectral shape. 
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Fig. 2. — EGB from SNe in star- forming galaxies. Results in the left panel assume pure 
luminosity evolution. Results in the right panel assume pure density evolution. The dashed 
blue line shows the contribution from CC SNe; the solid blue line shows the contribution 
from la SNe; and the black lin e plots the total con tribution from both CC and la SNe. The 
Fermi data are obtained from I Ab do et al.l (j2009al ). 



4.2. Pure Density Evolution 

For pure density evolution, only the galaxy density evolves with redshift while the galaxy 
luminosity does not. Therefore, the star-formation rate ip in a galaxy also remains constant, 
and the evolution in the cosmic star-formation rate will purely depend on the growth of 
the galaxy density. Hence, in the case of pure density evolution, Lho.z+Az = ^Ha,z- With 
similar calculations as those in the case of pure luminosity evolution (Appendix |XJ), one can 
find that (M gas S ) is also independent of the choice of the delay-time function. Additionally, 
(M gaS)S ) oc (1 + z)P in the case of pure density evolution. 

Results for the case of pure density evolution are shown in the right panel of Fig. HJ 
Again, the solid blue line and the dashed blue line represent the EGB from Type la and CC, 
respectively. The black line shows the combined gamma-ray contribution from both Type 
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Ia and CC events. Similar to the results of pure luminosity evolution, the EGB from Type 
la SNe is lower than that from CC events because of the lower Type la rate. Moreover, 
the predicted EGB emission is lower if we assume pure density evolution instead of pure 
luminosity evolution. As discussed in ([Fields et al.ll2010l ). this is because a typical galaxy's 
gamma-ray luminosity L 7 oc -?/>M gas ; in the pure luminosity evolution case, both factors are 
enhanced at early times, while in the pure density case this nonlinear boost is not present. 



5. The Extragalactic Gamma-ray Background from Type la Supernovae in 

Quiescent Galaxies 

Following a similar procedure to § HJ we will now discuss the EGB from cosmic rays 
accelerated by Type la SNe in quiescent galaxies. We again express the EGB luminosity 
density £ 7 a q (Eq. |9]) in the following form to describe the physics inputs from the average 
gas mass (M gaS) Q) and the cosmic la rate in quiescent galaxies 7Zu,q, 

/"la _ e Q ^Ia,Q pMW (p \ ( M \ 

7.Q ~~ c RMW , , pMW 1 if°^77^ em ' 1 ™ \ iW gas,Q/) y lU ) 



where 



/HST \ — (-^gas.Q -^Ia,Q ^galaxy,Q) /-, _\ 

{M gas> Q) = — r (17 J 

\-ttla,Q ^galaxy,Q/ 



dn 



,Q ' z ) dM 



dn 



dM*, 



and 7^i aj Q = (i?i a ,Q %aia xy ,Q)- Unlike the star-forming galaxies, where both (M gas S ) and 
7li a ,s can be related to the observable Ha luminosity, it is easier to connect both (M gaSj Q) 
and Hia,Q to the total stellar mass M^q in a quiescent galaxy. 



For the cosmic Type la rate in quiescent galaxies, we adopt the results of ISullivan et al. 



( 120061 ). which link Type la rates to M* q directly. These authors assume a bimodal delay-time 
distribution and decompose the la rate into two groups: the long-delay time and short-delay 
time. In their model, the short-delay time group simply traces the star-formation rate, while 
the long-delay time group has a constant probability for all delay times, i.e., A(r) = constant. 
Therefore, the Type la rate in a galaxy can be written as 

R la = A M* q + B?p. (19) 



The star formation rate ip ~ in a quiescent galaxy, thus -Ri a> Q = A M^q, where M* q is 
the total stellar mass created throughout the star formation history in the quiescent galaxy. 
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Sullivan et all (120061 ) estimated A = 5.1 x 10 14 yr 1 Mq 1 in quiescent galaxies based on 
measurements of the Type la rate in the Supernova Legacy Survey (SNLS). 



According to the observational results in lPannella et al.l ( 120091 ). the stellar-mass function 



dn 



of early-type galaxies evolves only slightly with redshift. Therefore we simply assume 

shown 



th e same stellar-mass fu nction throughout all redshift. Also, we find that the 
in 



dn 



dM*,Q 

Pannella et al.l ( 20091 ) can be roughly fitted by the following function □, 



dn 



d(log 10 M* iQ ) 



C m exp 



log* (M», Q //x) 
a 2 



(20) 



where C m = 2.05 x 10" 3 Mpc" 3 log" 1 M , fx = 10 ia7 M Q , and a m = 0.77. Equation 
(the stellar mass function for quiescent galaxies) and Eq. ( Tl9l) (Type la rate in a quiescent 
galaxy) give the cosmic la rate in quiescent galaxies 



K 



Ia,Q 



dM+ a R 



dn 



Ia,( 



1.53 x 10^ 5 yr 



dM^ tQ 
- 1 Mpc~ 3 . 



Here the constant rate follows from the nearly redshift-independent quiescent stellar mass 
function. 

Most of the gas content in quiescent gala xies appears to be i n the form of diffuse 
hot gas and can be observed in the X-ray (e.g., iForman et al.lll985t ICanizares et al.l 119871 : 
Bregman et al.lll992l ). However, large uncertainties exist in estimations of the mass of hot 
gas. Some studies suggest t hat most of the quiescent galaxies are gas-poor (e.g.. iDavid et al. 



2006; Fukazawa et al. 



20061 ), while other studies imp l y that there can be sig nificant amount 



of gas in these galaxies (e.g. iJiang fc Kochanekll2007l ; Humphrey et al.ll201ll ). Current theo- 
retical models suggest that the gas could extend out to a few hundred kpc from the center. 
The total gas mass of a single galaxy can increase by orders of magni tude if one includes 
the gas at large radii ( [Humphrey fc Buotd I201CH : Humphrey et al.l 1201 ll ) Therefore observa- 
tional measures of gas mass can vary even for a single galaxy, depending on whether or not 
the observations enclose a large enough radius to include gas at large distances. If hot gas 
extends to large radii, this could serve as a reservoir of targets for cosmic rays accelerated 



2 This fitting function is based on observations given in Fig. 7 in iPannella et al] ([20091 ) . which contains 
measurements of stellar-mass functions at different redshift bins and environments. We choose the data set in 
the lo west redshift bin of 0.25 < z < 0.55 in medium-dense environment (logiop 



Pannella et al 



—2.75; see 

(2009) for the definition of environmental density p and more details) to perform a \ 2 fitting. The reduced-%' 
of our fitting function is ~ 0.22. 
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by Type la SNe, provided the cosmic rays also reach large radii. This too is uncertain, and 
will depend on the cosmic-ray confinement versus escape properties of these galaxies. 

Because of the large uncertainties in the gas content in quiescent galaxies, we illus- 
trate a range of EGB predictions for these objects, based on three different gas amounts. 
Since quiescent galaxies are dominantly early-type, we will use estimates of the gas mass in 
ea rly-type galaxie s for t he amount of gas in quiescent galaxies. The stellar-mass function 



m 



Pannella et al.l (120091 ) (Eq. [2D]) gives the average stellar mass in early-type galaxies to 
be (M^q) = 5.17 x 10 11 M Q , which can be converted to the total gas mass by multiply- 
ing a gas-to-stellar mass ratio M gas /M*, i.e., (M gaSiQ ) = (M gaSiQ /M* iQ ) (M* iQ ). The three 
models we adopted for different gas amounts correspond to different gas-to-stellar mass ra- 
tios Mgag q/M^ Q for early- type galaxies. Table [T] summarizes the gas-to-stellar mass ratios 
^gas,Q,M^*,Q for the three models we adopted, as well as the corresponding (M gas Q). Gas 
Model 1 estima tes the gas-to-stella r mass ratio based on the stellar-mass fraction of the total 
halo mass from Jiang fc Kochanek ( 20071 ) §. This model gives the highe st gas-to-stell a r mas s 
ratio of all three models. Gas Model 2 and 3 are both adopted from iDavid et al.l ( 120061 ). 
in which the authors reported the gas-to-stellar mass ratio for more luminous (Gas Model 
2) and less luminous (Gas Model 3) early-type galaxies. Note that both of the gas mass 
and the la rate in a quiescent galaxy are constant with redshift, as a result of assuming 
a non-evolving stellar -mass function, consistent with the early-type galaxy observations in 



Pannella et al. 



(I2009J). 



Table 1: Summary of different gas models adopted in the EGB calculations. 



Gas Model 


M gas , Q /M*, Q 


( M gas,Q) 


Reference 


1 


2.14 


1.11 x 10 12 M 


Jiang & Kochanek (2007) 


2 


0.01 


5.17 x 10 9 M 


David et al. (2006) 


3 


0.001 


5.17 x 10 s M 


David et al. (2006) 



3 Ijiang fc Kochanekl (|2007l ) found that the average stellar-mass fraction of the total halo mass in early- 
type galaxies is M+/Mtot ~ 0.026 or 0.056 based on different assumptions of the halo mass dynamics. 
Both of these numbers are significantly lo wer than the cosmological baryon-to-mass ratio fli,/n m ~ 0.176 
measured by WMAP (jSpergel et al.l 120071) . If we assume that the baryon-to-mass rat io in a galaxy can 
be we ll represented by the cosmological ratio, i.e., Mb/M to t ~ ^b/^m, the result from Ijiang fc Kochanek 
(|20071 ) implies a large amount of gas mass in early-type galaxies, which can be estimated by M gas q = 
(M baryoniQ - M*, Q ) ~ M*, Q ( M*fuL ~ ^ The values of M */ A/ tot ~ 0.026 and 0.056 correspond to 
M gaSy Q = 5.77 M^q and M gaS! Q = 2.14 -M^q, respectively. Here we adopt the latter number to be more 
conservative in our estimation. 
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The red curves in Fig.[3]plot the EGB estimation from Type la SNe in quiescent galaxies, 
with different line styles correspond to estimations from different gas amounts (solid line: 
Gas Model 1, dashed line: Gas Model 2, dotted line: Gas Model 3). The EGB emissions 
from SNe in star-forming galaxies are plotted as blue curves for comparison. The black 
curves in Fig. H] plot the total EGB emissions from SNe in both star-forming and quiescent 
galaxies for different gas models. For both Fig. [3] and HJ the left panel plots results under 
the assumption of pure luminosity evolution for the star-forming galaxies. The right panel 
shows the EGB predictions assuming pure density evolution for the quiescent galaxies. 

The estimated EGBs shown in Fig. [3] are linearly proportional to the average galaxy gas 
mass (M gaSj Q) and thus to the adopted gas-to-stellar mass ratio, as required by Eq. (fT6l) . Gas 
Model 1 predicts the highest contribution in EGB from quiescent galaxies, where the result 
is significantly higher than the EGB emission from Type la events in star-forming galaxies, 
while Gas Model 2 and Gas Model 3 suggest much lower EGB emission from quiescent 
galaxies. In our estimation, the two important factors that affect the overall EGB are the 
average gas mass and the cosmic SN rate. In general, the cosmic Type la rate is about 
a factor of 5 — 10 smaller than the cosmic CC SN rate. Additionally, the cosmic Type la 
rate in all quiescent galaxies averaging over the entire redshift range is lower than that in 
star-forming galaxies by around a factor of three. 

As for the gas mass in quiescent galaxies, Gas Model 1 assumes the gas amount to be 
about two orders of magnitude higher than that in star-forming galaxies, while Gas Model 
2 and 3 assume gas amount to be similar or one order of magnitude lower than that in star- 
forming galaxies, respectively. Combining these two factors (gas mass and Type la rate), we 
would expect the EGB from quiescent galaxies to be about 30 times larger (for Gas Model 
1), 3 times smaller (for Gas Model 2), or 30 times smaller (for Gas Model 3) than that from 
Type la in star-forming galaxies. Because of the EGB estimation is very sensitive to the 
gas amount, measurement of the EGB could put constraints on the gas mass in quiescent 
galaxies. 

6. The Uncertainties in the Extragalactic Gamma-ray Background Analysis 

6.1. Star- forming Galaxies 

The gamma-ray emission from normal star-forming galaxies is the most constrained 
among the three galaxy classes that are considered possible SN-induced EGB sources (i.e., 
starburst galaxies, normal star- forming galaxies, and quiescent galaxies). The main uncer- 
tainties in the EGB prediction of star-forming galaxies come from four factors, as described in 
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Fig. 3. — EGB from SNe in both star-forming and quiescent galaxies. Results in the left 
panel assume pure luminosity evolution for star-forming galaxies. Results in the right panel 
assume pure density evolution for star-forming galaxies. The dashed blue line shows the 
contribution from CC SNe; the solid blue line shows the contribution from la SNe in star- 
forming galaxies; the red lines show the contribution from la SNe in quiescent galaxies based 
on different gas models (solid red: Gas Model 1, dashed r e d: Gas Model 2, dotted red: Gas 
Model 3). The Fermi data are obtained from I Ab do et al.l ( 2009a ). 



Fiel ds et al. 1) uncertainty in the pionic gamma-ray production rate T^^^E, 



which is ~ 30% (lAbdo et al. 



2009 cJ), (2) uncertainty in th e normalization of the Galactic SN 



which is ~ 40% (IRobitaille fc Whitney! 120101) . (3) uncertainty in the luminosity 



rate R^^^ ^ 

scaling in (M gas S ), which is ~ 25% (jFields et a l 
ization of the cosmic SN rate 7£sn,s> which is 
cosmic CC SN rate TZ CC ~ (1.0 ± 0.2 



201ol ). and (4) uncertainty in the normal- 



16% resultin g from the uncertaint ies in the 



x 10~ 4 yr" 1 Mpc~ 3 



Horiuchi et al. 



cosmic la rate K Ia ~ (0.25 ± 0.05) x 10" 4 yr" 1 Mpc~ 3 iHoriuchi fc Beacomllioiol ). The total 
uncertainty in the EGB prediction will then be ~ 10 ±0,25 . 



2009) and the 



The upcoming large synoptic surveys, such as the LSST, will provide novel information 
in both the cosmic SN rate and how they depend on the galaxy types out to high redshift. 
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Fig. 4. — Total EGB from la and CC SNe in both star-forming and quiescent galaxies. 
Results in the left panel assume pure luminosity evolution. Results in the right panel assume 
pure density evolution. The Fermi data are obtained from I Ab do et al.l (j2009al ). 



Within one year of observation, LSST is expected to detect ~ 10 5 SNe out to z ~ 1 and thus 



achieve a statistical prec ision of less than a few percent in the cosmic SN rate (jLien fc Fields 
20091 ; iBailey et al.l 120091 ) . Hence, LSST will essentially remove the uncertainty from the 
cosmic SN rate in the EGB analysis and make the EGB a better tool for studying cosmic 
rays and gamma-ray physics. Moreover, such a large SN population will provide excellent 
statistics for study of how SN rates evolve as a function of SN type (and sub- type), host 
galaxy type and star- formation rate, and cosmic environment. Thus it will be possible to 
directly measure the Type la rate in quiescent galaxies and its evolution with redshift. 



6.2. Quiescent Galaxies: Hot Gas and Cosmic-Ray Propagation 

For quiescent galaxies, many characteristics related to their gamma-ray emissions are 
poorly understood. The Type la SN rate in these systems could be uncertain up to a factor 
~ 2, particularly at z > 1. We have also seen that the gas mass in these systems is even more 
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poorly known. Published estimates of the gas content in quiescent galaxies vary by orders of 
magnitude. Observation s that use X-rays a s tracers of hot gas fo und significantly less gas in 



quiescent galaxies (e.g., iDavid et al.ll2006t iFukazawa et al.l 120061 ). However, several studies 



imply a much higher gas content in quiescent galaxies than previo usly thought, when using 
gravitational lensing and modeling that in volved dark matter (e.g.. |Jiang fc Kochanekl 120071 ; 
Capelo et al.ll2010t Humphrey et al.ll201ll ). It is difficult to determine what are the causes of 
the large discrepancy of gas quantity from different analyses, because these studies adopted 
different observational and modeling techniques with different galaxy sample. One possible 
reason for the large variation in gas content measurements might come from w hether or not 



one includes gas at larger radii, which are predicted by s ome theoretical models (jCapelo et al. 
20 id iHumphrev & Buotdlioia Irlumphrev et alJboilh . 



If gas content does extend to larger radii in quiescent galaxies, it is crucial to understand 
the cosmic-ray propagation in quiescent galaxies to determine the confinement volume of 
cosmic rays and how likely they can interact with gas at larger radii. Unfortunately, cosmic- 
ray propagation remains poorly understood. Until now, most of the studies have been 
focused on the Milky Way or spiral galaxies and not so much on quiescent galaxies. For spiral 
galaxies, both observations and th eoretical modeling suggest that most of the cosm ic rays are 
confined within ~ kpc of the disk (jStecker fc Jonedl 19771; IStrong et al.ll2000U2004 ). However, 



this might not be the case for quiescent galaxies. iHein &: Spanierl (120081 ) simulate cosmic- 
ray propagation in elliptical galaxies with the presence of diffusion, as well as adiabatic 
losses. They apply their formalism not to SN sources but rather to acceleration due to 
an M87-like relativistic jet; thus their detailed calculations are not applicable to our case. 
Nevertheless, their general finding is that cosmic rays expand into a larger volume in an 
elliptical galaxy than in a spiral galaxy, due to the larger minor axis in an elliptical galaxy. 
Additionally, they argue that adiabatic losses are much more i mportant in elliptical galaxies 
and indeed dominate over escape losses. IHein fc Spanierl (120081 ) thus conclude that it is likely 
for elliptical galaxies to be extended gamma-ray sources. 

There are two factors that are important for determining the probability of interactions 
between cosmic rays and gas, especially at large radii where the density of gas and cosmic 
rays are likely to be smaller: (1) the mean free time for pion production for each cosmic-ray 
particle r. 



pp— >7T 



(I 



^ 1 Gyr (0.1 cm 2 /n gas ), where (r pp ^. n o is the cross section 
of pion production and n gas is the number density of gas particles, and (2) time r esc that 
a cosmic ray particle takes to propagate through the galaxy before escape. In the limit of 
cosmic rays travel through a (non-magnetized) galaxy radially, that is, without any diffusion, 
the escape time is very short and probability of interactions between cosmic rays and gas 
particles is small. On the other hand, in the presence of disordered magnetic fields extending 
to large radii, cosmic ray propagation will not be radial but diffusive instead, which will make 
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the escape time much larger and increase the chance of cosmic-ray interactions with ISM. 



Clearly there are fundamental uncertainties (and opportunities) in a realistic treatment 
of elliptical galaxy cosmic rays, but even in our simplistic picture the parameters A esc and 
€q are not well-constrained. Although we treated these two quantities to be the same for 
both Type la and CC SNe in all environments due to limited knowledge, it i s possible, and 
even l ikely, that these numbers are different in quiescent galaxies. In fact, iDorfi fc Voelk 
( 119961 ) have suggested that th e emciency in quiesce nt galaxies is at least 10 times lower than 
that in star-forming galaxies ( Dorfi fc Voelk! Il996[ ). which coul d lower our prediction of the 
EGB from quiescent galaxies by a factor of 10 or even larger. iTang fc Wangl ( 120051 ) model 
the evolution of SN remnants in low-density hot media and reach similar conclusion that the 
SN heating in such an environment is subtle and cover a large region because of small Mach 
numbers. Likewise, a smaller escape path length, i.e., a weaker cosmic-ray confinement, can 
also decrease our EGB estimation in quiescent galaxies. Note that adopting different values 
of e and A esc would change the cosmic-ray spectrum, and hence also change the corresponding 
gamma-ray spectrum. 

To summarize, current studies imply that cosmic-ray acceleration might be more diffi- 
cult in quiescent galaxies than in star -forming galaxies du e to weaker shocks in low-density 
and high-temperature environments (jDorfi fc Voelk! Il996l ; ITang &: Wangl 120051 ) . However, 
quiescent galaxies might have larger confinement v olume than spiral galaxies because they 
have larger semi-minor radii (jHein fc Spanierll2008l ). Further observational and theoretical 
study of the gas and cosmic-ray content of quiescent galaxies clearly are needed in order to 
pin down the gamma-ray production of these objects even to within an order of magnitude. 



Conclusions 



We have calculated the EGB contribu tion from Ty p e la S Ne in both star-forming and 
quiescent galaxies, extending the work of Fields et all ( 120101 ) . For star- forming galaxies, 
most of the gamma-ray emission comes from cosmic rays accelerated by CC SNe. This is 
mainly because there are about five times more CC SNe than Type la events in star-forming 
galaxie s. We find t h at th e net EGB contribution from both SN types is almost the same 
as the iFields et all (120101 ) predictions that only included CC events. Our model allows 
for addition of cosmic Type la explosions, but also includes these in the cosmic-ray/star- 
formation ratio, which we normalize to the Milky Way values. Both factors change by nearly 
the same amount, so that the addition of Type la events causes almost no net change in the 
EGB prediction. 
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We also point out that Type la events in quiescent galaxies make a unique contribution 
to the EGB, because these systems lack CC events. We show that the EGB from Type 
la events in quiescent galaxies is highly sensitive to the gas amount in quiescent galaxies, 
which is still poorly known. Based on different gas models adopted, the EGB from Type 
la SNe in quiescent galaxies can vary from two orders of magnitudes higher to an order of 
magnitude lower than those produced by Type la SNe in star-forming galaxies. Therefore 
quiescent gal axies can be dominant sour ce of the EGB if there e xist a large amount of gas, as 



suggested by IJiang fc Kochanek! (120071 ) and Humphrey et al.l (120111 ) . The measurement of 



the EGB will provide useful constraint on the gas amount in quiescent galaxies. Additional 
uncertainties in the cosmic-ray acceleration efficiency and confinement could also change the 
EGB emission in these systems by a few orders of magnitude. Hence the EGB can also 
provide limits on the these two quantities in quiescent galaxies. 

It is thus important to understand the characteristics of cosmic SNe of all types, in order 
to correctly predict their contribution to the EGB. We conclude that the large SN sample 
provided by LSST will offer critical information about the cosmic SN rate for both CC and 
la events, and their dependence on galaxy types out to high redshift. 

The Fermi detection of the EGB contains crucial information about the extragalactic 
gamma-ray source spectrum. Particularly, it can provide an important probe to the gas 
amount, the cosmic-ray acceleration efficiency and the cosmic-ray confinement in quiescent 
galaxies. With our knowledge about SNe increasing rapidly as future synoptic surveys come 
online, the EGB contribution from SNe in galaxies can possibly be disentangled from other 
source candidates. 
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A. The Delay-Time Distribution of Type la Supernovae and Detailed 

Calculation of (M gas s) in Equation [TOl 



The delay time of each Type la SN can differ from ~ 0.1 Gyr to ~ 10 Gyr (IMannucci et al. 



2005 



2006 



Maoz et al. 



Scannapieco fc Bildsterfe ood: bal-Yam fc Maojl2004 ; IMannucci et al. 2006; Sullivan et al 



20111 ). Observationally, the delay times of Type la SNe are usually stud- 
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ied via comparison between measurements of the cosmic Type la SN rate and the cosmic 
star-formation rate, and are usually described by some functions of delay-time distribution, 
which describes the probability of a Type la event with a specific dela y time. Current pro 



posed delay-time distri butions have inclu ded a single power law (e.g., iHoriuchi &: Beacom 



2010l : iGraur et al.lhoilh , a Gaussian (egJstroker et al.l l2004t bahlen et alll2008h , and a bi 



modal distribution (e.g., ISullivan et al.l 120061 ) . The difficulty in determining the delay-time 
distribution mainly comes from the uncertainty in the cosmic la SN rate measurements. 
Fortunately, the value (M gas S ) is independent of the choice of delay-time distribution, as we 
show in the following derivation. 



Based on the relation between the gas mass M gaS) s and the star-formation rate ip (Eq. fTB"]) . 
and also the connection between ip and Lh«, z , i>e., ip(Ln a ,z)/(l M Q yr _1 ) = L Haz /(1.26 x 
10 34 W), the M gas> s for a galaxy at redshift z can be directly linked to the observable H a 
luminosity L Hq , z by M gas>s = 2.8 x 10 9 M (1 + z)'? { 12 ^Ia w )". Additionally, the Type 
la rate can also be related to Lh q , z by 



#ia,s oc / L Ha (t-r) A(r) dr, 



(Al) 



where L Ha {t — r) = L Ha z+ A z (t, r), which is the H a luminosity measured at some earlier time 
t — r or larger redshift z + Az. Therefore (M gaSj g) can be expressed in terms of Lh q , z , 



(M gaS)S ) CX 



/ dL Ha:Z (l + z)~P (L Ha , z ) u (/J L Ha , z+ Az(t,r) A(r) dr) 



dn 



galaxy. S 



JdL RatZ (Jj L HQ , Z+Az (t,r) A(r) dr) 



(A2) 



This equation expresses only the redshift-dependent terms and we will further calculate 
how (M gas S ) evolves with redshift in the case of pure l umino sity evolution and pure density 
evolution, respectively, in IA.1I and IA.2I iFields et al.l ( 120101 ) shows that the local value of 

(M gas , s ),=o = 6.8 x 10 9 M . 



A.l. Pure Luminosity Evolution 



In the case of pure luminosity evolution, the star-formation rate in each galaxy traces 
the general evolution of the cosmic star-formation rate, i.e., = p* (£+Ag) Therefore 

from the directly proportional relation between the star-formation rate if) in a galaxy and the 
galaxy H a luminosity I/h q , z , one can trace the evolution of the H a luminosity via the history 



of cosmic star-formation r ate, which is well-known out to redshift z ~ 1 
Hopkins fc Beacom 2006 . and references therein). That is, L "" A ' 



Hopkinsll2004 



p+jz+Az) 
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Therefore the galaxy luminosity at different redshifts can be found by 



p*(z + Az) _ pS~ t ) 

^Ha,z+Az — ^Ha.z . , \ = ^Ha.z . , , N • 

With this relation, (M gas S ) in Eq. (IA2[) can be simplified to 



(A3) 



(M gaStS ) oc 



P*(t-r) 

Mt) 



A(r) dr) 



dn 



galaxy ,S 



dL 



Hq.z 



till. 



galaxy, S 



galaxy ,S 
^Ha.z 



/ dLu a ,z (Lnc 



W+l rf "galaxy,S 



Ha,z 



dn 



galaxy, S 



(A4) 
(A5) 
(A6) 



Therefore (M gaS) s) is independent of the choice of delay-time distribution. 



The assumption of pure luminosity evolution implies that ra* i2 in the Schechter function 
remains constant and L* iZ in the Schechter function evolves as p\ (Eq. IA3I) . Hence the 
redshift dependence of (M gaS) s) can be further calculated using the Schechter function, 



(M gas , s ) ex (1 + z)-* L\ 



^Ha,z \u)+l 



:J-a g— Lho,,z/L* 



L*, z V i 
'Lh q .z\ /Lh 



g~^Ha,z/^* 



(A7) 



*,z=0 



P*(^ = 0) 7 j L ^d(^ 



iHa.z 71* 



oc (1 + z) 



(A8) 
(A9) 



where L max is the maximum lumino sity for star- f ormin g galaxies, which corresponds to the 
maximum star formation defined in iFields et al.l (120101 ) . Galaxies with luminosities greater 
than L max are considered starburst galaxies and are not included in this calculatio n . Ad- 
ditionally, we adopt the cosmic star-formation rate described in iHoriuchi et al.l (120091 ) 
based on current observations. Note that because the factors related to delay-ti me distribu- 



tion c anceled out, this result turns out to be the same as the one obtained in IFields et al. 

(boioh . 
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A. 2. Pure Density Evolution 

In the case of pure density evolution Lh«, z +a z = -^Ha.z as discussed in § 14.21 Thus in the 
Schechter function, L* )2 remains constant while n* iZ evolves as p*. With similar calculations 
shown in the case of pure luminosity evolution (IA.1I) and adopting the Schechter function for 
^gaxy.s ^ we can derive rec lshift evolution of (M gaSi g) in the case of pure density evolution: 



/ dL HaiZ (l + z)-P (L Ha , z r (Jo L na>z A(r) dr) 



dn 



galaxy,S 

dLuc 



< J 'W> oc - 7., . . „„ lTO , s (MO) 



/ cLL H «, z (Jo L HoijZ A(r) dr) 



dL 



Ha.z 



l + z)-/ 3 Cg A(r) dr) JcjL HQ , z 

(f Q A(r) rfr) /rfL Ha , z L Ha , z %^f 



dn. 



galaxy, S 



(All) 



(A12) 



> z d,L HctyZ 



oa{l + z)- p . (A13) 
Again, because (M gas g) is indepe ndent of the choice of delay-time distribution, the result is 



identical to the one calculated in iFields et al.l ( 120101 ) 
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